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INTERPOLATION AND SMOOTHING 
OF EXPERIMENTAL DATA 
WITH SLIDING POLYNOMIALS 

By W. M. Snyder' 

ABSTRACT 


This paper presents a comprehensive overview of the analytical structure 
of sliding polynomials. Emphasis is placed on explanation so that researchers 
with limited background in numerical analysis can apply the methods. 
Presentation is generalized to assist readers from various disciplines in 
preparing data for interpretation. Polynomials of increasing complexity were 
previously presented in separate publications. The intent here is to develop 
an appreciation of the common elements underlying three levels of complexi- 
ty, starting with simple linear forms. Continuity of the first derivative of 
interpolation arcs is added first, followed by extension to continuity of the 


second derivative. The three forms 
operators using least squares. 

INTRODUCTION 

Interpolation and smoothing are routinely re- 
quired to reduce and process experimental data. 
Choice of a particular method depends on the 
amount of data and on awareness of available 
methods. DuChateau et al. (1972) and Kimball 
(1974) dealt with interpolation and smoothing 
techniques that require substantial amounts of 
experimental data. Such data may be generated 
either in an extensive cyclic time domain or by 
repetition in a limited time or space domain. 
Savitzky and Golay (1964) discussed smoothing 
and differentiation of data using convoluting 
functions in weighted moving average tech- 
niques. Again, availability of substantial a- 
mounts of data is implied, generated by basically 
continuous physical experimentation. The author 
(1961, 1962, 1967a) showed how sliding polyno- 
mials could be used for interpolation and smooth- 
ing. This technique can be used with more 
limited amounts of experimental data. 

' Hydraulic engineer, Southeast Watershed Research 
Laboratory, Agricultural Research Service, U.S. Depart- 
ment of Agriculture, Athens, Ga. 30601 . 


are then transformed to smoothing 


This publication draws together and expands 
the author’s previous work on sliding polyno- 
mials. The construction of sliding polynomials is 
reviewed and analyzed. Common elements and 
concepts can be found in a hierarchy of polyno- 
mials. Tracing these common elements through 
three levels of complexity will be useful to 
researchers in many disciplines who must pro- 
cess observational or experimental data. Realiza- 
tion of the full potential of the sliding polynomials 
will aid in choosing interpolation and smoothing 
techniques best suited to specific data processing 
requirements. 

Understanding the hierarchial structure of the 
polynomials can begin with review of simple 
linear interpolation. Extension of this simple 
form to piecewise linear functions establishes a 
basis for understanding more complex forms. 
Addition of a requirement for mathematical 
continuity of a first derivative is a logical next 
step. Further extension follows with require- 
ment for mathematical continuity of a second 
derivative. All of these three piecewise forms are 
utilitarian. While a basic knowledge of mathe- 
matical and numerical analysis is assumed, em- 
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phasis is placed on explanation to aid the reader 
in selection of the form suited to his purposes, 

INTERPOLATION 
WITH POLYNOMIALS 

Interpolation is the estimation of values be- 
tween given discrete values of a property. Some 
function must be assumed. This function must 
give values of the property dependent upon 
values and locations of the given discrete points. 
A one-to-one correspondence exists between 
mathematical degrees of freedom in the function 
and the number of given discrete points. Data 
points through which the function is required to 
pass will be called base points. 

In the following discussion, basic polynomial 
properties and the functional forms of sliding 
polynomials are explained. Then, the use of these 
concepts in interpolation is described. 

Basic Polynomial Properties 

Polynomials of any degree could be assumed 
for interpolation. Beyond the fifth power, how- 
ever, the algebraic manipulations become tedious 
and error prone, A fifth -power polynomial, 
evaluated from six or more coordinate points, is 
sufficient to produce smoothly continuous rates 
of change through experimental data. Therefore, 
this presentation does not go beyond fifth-power 
polynomials. 

The fifth-power polynomial has the conven- 
tional form 

y—a+bx+cx^+dz^-hex^-hfic^. ( 1 ) 

Coefficients a through / represent six mathe- 
matical degrees of freedom which can be satisfied 
by reference to six given discrete points. These 
concepts are developed under “Six-Point Inter- 
polation." 

The first derivative of equation 1, the rate of 
change of y, is 

dy/d3:^b^2cx-^^dx'^'^iex'^’^b/x\ ( 2 ) 

Higher order derivatives are given in equa- 
tions 3 through 6. 


d^y/dx^^2c-\-Qdx-\‘l2ex 

(3( 

d^y/dx 24 ea: -f 60/a; 

(41 

d'^y/dx'^—2Ae-hl20fx. 

(5) 

d^y/dx^^l20/. 

(6) 


Note the fundamental continuity of form in 
equations 1 through 5. Any one of these can be 
used as an interpolation function. Properties of 
the function for estimating rates of change are 
given by derivatives of order higher than the 
chosen functional form. Properties of the func- 
tion for integration are given by derivatives of 
order lower than the chosen functional form. For 
example, equation 5, the straight line, has a 
constant rate of change given by equation 6 and a 
parabolic rate of accumulation given by equation 
4. In the discussions of two-point, four-point, and 
six-point interpolation, equations 5, 3, and 1 are 
used to develop special polynomials. They act as 
flexible splines that pass through extended num- 
bers of coordinate points. Since symmetry is 
required, only polynomials with even degrees of 
freedom can be used in this method. 


Sliding Polynomials 

Different terminology may be used to describe 
the mathematical techniques of passing a contin- 
uously transhaping function from left to right 
through some set of discrete values, or base 
points. As the function moves to the right, base 
points must drop out on the left side of the 
function span, and new points must be picked up 
on the right side of the function span. This 
concept can be called continuous transforming, 
flexible spline, or sliding polynomial. The term 
“sliding polynomial" is considered most suitable 
in the development that follows. 

First, consider the schematics in figure 1. They 
are sufficient for a simple situation where inter- 
polation is always within one interval between 
base points. 

Note the set of coordinate points YBO through 
YBQ in figure la. These coordinate points, spaced 
at uniform intervals, wUl be used as base points 
through which interpolation polynomials are re- 
quired to pass. The designation PBO, for exam- 
ple, means a base point in y at a value of x equal 
to zero. Uniform spacing of the base points is not 
absolutely necessary, but it does greatly simplify 
the algebraic forms. 

Now observe in detail the interval between 
YB2 and FBS. Figure lb is the schematic 
representation for passing a straight line be- 
tween the two points. Interpolation can be made 
on this straight line. Figure Ic shows a cubic 
equation passing through YBl, YB2, YB3, and 
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0 12 3 4 5 6 


a. DEFINITION OF BASE POINTS 

Y02 YB3 

O— — O 

y= 0 + bx 

b. TWO- POINT INTERPOLATION 


V I 

y=a + bx + cx2 + dx^ 
c. FOUR- POINT INTERPOLATION 


y=a + bx+ cx2 + dx3+ex‘'+fx® 
d. SIX- POINT INTERPOLATION 

FIGURE 1 . — Definition of Interpolation systems. 


on the left and picking up FB5 on the right. The 
phasing out and phasing in of the base points is 
gradual, instead of abrupt, as in the linear 
segments. 

Figure 2c is an extension of concepts in figure 
2b. The higher order polynomial phases out YBO 
and phases in 756 as the a;-point leaves the 
YB2-YBZ arc and passes through x=3 into the 
753-754 arc. Special forms for this method are 
shown under "Six-Point Interpolation." 

Two-Point interpolation 

Two-point interpolation is the familiar and 
conventional linear interpolation. The method is 
intuitively understood and used by many per- 
sons. This intuitive understanding is used here as 
the basis for presentation of more complex 
polynominal methods. To prepare for this pres- 
entation, the familiar linear interpolation is cast 
into algebraic forms which will be needed later. 
Only two base points at a time are used in linear 
interpolation. In figure la, the two points are 
752 and 753. The two parameters in the 
associated equation can be evaluated in terms of 
the base points. The matching of number of 
parameters and number of base points should be 
noted. 


754. Interpolation is to be limited to the interval 
between 752, and 753, the same interval cov- 
ered by the straight-line example. Figure Id 
shows a fifth -power equation passing through six 
points. Interpolation again is to be limited to the 
interval between 752 and 753. 

A more realistic, and more complicated, prob- 
lem is interpolating in adjacent intervals (fig. 2). 
Figure 2a shows interpolation on two adjacent 
linear segments. Visualize a point of interpola- 
tion on the x-scale moving from 752 to 753. 
Values of y can be read from the straight line 
connecting 752 and 753. When the point passes 
through 753, a transformation is necessary. 
Point 752 is dropped on the left, and point 754 is 
added on the right. As the x-point moves through 
x=3, interpolation changes abruptly from the 
line through 752 and 753 to the line through 
753 and 754. 

A smoother transition through 753, utilizing 
an arc through 752 and 753 and one through 
753 and 754, is shown in figure 2b. The 
technique to suppress the angle between linear 
segments at 753 is developed under “Pour-Point 
Interpolation." It is based on dropping out 751 



a. CONNECTED LINEAR SEGMENTS 



b. SLIDING PARABOLA 



c. SLIDING FOURTH POWER POLYNOMIAL 

FIGURE 2. — Definitions for connected arcs. 
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Define a new scale variable z, with «=0 at 
x~2, and z=l atx=3. An equation to define the 
assumed linear interpolation function between 
YB2 and 753 can now be written: 

y=a+bz. <7) 

At i=0, y=YB2, and z=l, y=YB3. Substi- 
tuting these coordinate values of the two base 
points in equation 7 yields the simultaneous 
equations 

YB2=(i (8) 

YB3=a+b. (9) 

From equations 8 and 9, a=YB2, and b~YB3 
~YB2. Equation 7 may now be written in terms of 
the base points as 

y=YB2+(YB3-YB2)z={l-z)YB2+z-YB3. (10) 

The solution for equations 8 and 9 may also be 
written in the form of the matrix-vector equation 


1 

0 


yB2 


a 

-1 

1 


YB3 


h 


So-called interpolation coefficients can be de- 
veloped from special values of 752 and 753. 
First, set 752 equal to 1 and 753 equal to zero. 
Expansion of equation 11 gives a=l and 6=-l. 
The special form of equation 10 results in 

Vo=\~z. (12) 

Since 752 was set to unity, yu can be called the 
unit contribution of a base point at 2 = 0 , the left 
one of two base points. Next, set 752 equal to 
zero and 753 equal to 1. Expansion of equation 
11 gives a=0 and 6=1. The special form of 
equation 10 now is 

Vh-z. (13) 

Now, t/j can be called the unit contribution of a 
base point at 2 = 1 , the right one of the two base 
points. 

Equations 12 and 13 can be evaluated for 
values of z as in table 1. The values of y„ and 
are normally designated as CO and Cl, meaning 
the interpolating coefficient to be multiplied by 
the base value at 2=0 and by the base value at 
2 = 1 . The same meaning is expressed on the 
right in equation 10. For example, two-point 
interpolation at a spot one-fourth of the x-dis- 
tance from 752 to 753 is written directly, with 
coefficients from table 1, as 

y=(0.75)yB2-)- (0.25) KBS. ( 14 ) 

Equation 14, of course, is the familiar and 


conventional form for linear interpolation. 

If the two adjacent linear segments in figure 2a 
are considered, equations can be written directly 
for interpolating y , on the left segment and y^ on 
the right segment; 

y =YB2+iYB3-YB2)z, (15) 

and y =YB3-\-{.YB\~YB3)z. (16) 

With 2=1, equation 15 gives y, ~YB3. Also, 
with 2=0, equation 16 gives y., — YB3. This 
simply confirms that both segments pass through 
753, and the interpolation is continuous in the 
sense that there is no gap at the base point 753. 
However, dyi /(f 2 = 753 -752 and dy-i/dz=YBi 
-753 show that the slopes of the segments are 
different, constant, and not continuous. 


Table 1 — Two-point interpolation coefficients 


z 

Vn 

Vb 

0,10 

0.90 

0.10 

.25 

.75 

.26 

.50 

.50 

.50 

.75 

.25 

.76 

.90 

.10 

.90 


Four-Point Interpolation 

In geometrical concept, two-point interpola- 
tion can be visualized as a point being linearly 
displaced in y-scale as it moves along the 2-scalo 
from one base point to the next. This two- 
directional displacement generates a locus of 
interpolated values which lie along a straight 
line. Note particularly the concept of a point in 
continuous displacement generating a straight- 
line interpolation locus, and note that the line is 
one mathematical order up from the point. 

In four- point interpolation, the geometrical 
concept is a function that will be linearly tran- 
shaped in y-scale as a point of interpolation 
moves along the 2 -scale. This function is up two 
orders from the point, and is, therefore, a 
parabola. The locus of interpolated values is one 
•order higher than the function being transhaped 
and is a cubic polynomial. 

The geometrical concept on which four-point 
interpolation is based is further developed as 
follows: Visualize a parabola that passes through 
base points 751, 752, and 753 in figure 1. Next, 
visualize the parabola through 752, 753, and 
754. As the point of interpolation moves along 


4 



the x-scale from YB2 to F53, the interpolated 
value is to be read from a flexible parabola that 
transhapes linearly from the left parabola to the 
right parabola. 

The analytical concept of four-point interpola- 
tion requires specification of four mathematical 
constraints on which solution for the four coeffi- 
cients of the cubic polynomial, equation 3, are to 
be based. The mathematical specifications are 
based on commonality. The cubic equation is 
required to be common to the parabolas at the 
ends of the interpolation arc. This means simply 
that both the parabolas and the cubic pass 
through YB2 and KB3 (fig. Ic). In addition, the 
parabolas and the cubic are required to have the 
same slope, or the same first derivative, at YB2 
and F53. 

Next, specify a scale variate z, which takes 
values 0, 1, 2, and 3 from left to right in any 
four-point set of base points. Define the left 
parabola by 

\y=u-\-vz+wz'^ (17) 

and d\y/dz=v+2wz. (18) 

Solve for coefficients w, v, and w by passing 
through the three points «=0, iy~YB\\ »=1, 
,j/=FB2; and 2=2, ^y—YEZ, giving the values 


u =KB1, (19) 

v=-(3/2)FJ?H-2-yB2-(l/2)KB3. (20) 

and w=(l/2)KBl-yB2+(l/2)FB3. (21) 

Next define the right parabola by 

2 V—u-{-vz-\-u>z (22) 

and d2y/dz~v\-2wz. (23) 


Solve for new coefficients u, u, and w by 
passing the parabola through points z~l, 
2y=YB2,\ 2=2, zy—YB^-, and 2=3, 2y=FB4. 

The solution is 


«=3-yB2-3-yB3-|-FS4, (24) 

v=:-(6/2)ra2-|-4F53-(3/2)yB4, (25) 

and w=(l/2)YB2-YB3+{l/2)YBA. (26) 

Define the cubic equation by 

y—a+bz+cz +dz^ (27) 

and dy/dz^h-\-2cz'^i^z‘^. (28) 


The four mathematical constraints expressing 
commonality ^tween parabolas and cubic can 
now be stated as four simultaneous equations 
(29-32), Equations 17 hnd 18, with coefficients in 
equations 19, 20, and 21, are set equal to 


equations 27 and 28 for 2=1. Equations 22 and 
23, with coefficients in equations 24, 25, and 26, 
are set equal to equations 27 and 28 for 2=2. The 
resulting four simultaneous equations are 


m=a-l-6-fc-l-d, (29) 

(l/2)yB3-(l/2)rai=i>-)-2c-f3d, (30) 

753=«-)-26-(-4c-|-8d, (31) 

and (l/2)yfi4-(l/2)m=b-f4c-l-12d. (32) 


The solution to these equations is 

a=2.yBl-3-yB2-f3-FB3~ra4, (33) 

b=-i-YBl+m/2)YB2-8-YB3+{5/2)YBi, (34) 

<:=(5/2)m-7-FB2+(13/2)KB3-2-yB4, (35) 


and d=-{\mYB\-\-[3/2)YB2~(3/2)YB3MV2)YBL (36) 

This solution can be expressed concisely in the 
matrix-vector equation 


4 - 6 6-2 


YBl 


a 

- 8 19 -16 6 


YB2 


b 

6 -14 13-4 


YB3 

— 

c 

- 1 3 - 3 1 


YBi 


d 


Next, consider the property of the locus of 
interpolated values at the common point of two 
arcs. Specifically, in figure 2b, arcs F52-FS3 and 
YBZ-YBA obviously join at YBZ. Now at 753, 
interpolation arc 752-753 ends on the parabola 
through 752, 753, and 754, However, interpo- 
lation arc 753-754 takes off from this same 
parabola. The four mathematical specifications 
require the interpolation cubics to be common 
and tangent to the parabolas at the ends of the 
arcs. The cubics of the 7B1-TB4 set and the 
752-755 set are, therefore, both common and 
tangent to the same parabola at 753. The cubics 
of adjacent interpolation arcs are thus common 
and tangent to each other. 

Thus, it has been demonstrated that four-point 
interpolation based on a sliding, transhaping 
parabola generates values on a cubic polynomial. 
Adjacent cubic arcs are continuous and have 
continuous first derivatives as the interpolation 
locus passes through the base points. 

Figure 3 is an example of four-point interpola- 
tion. (Computer program A-1, appendix A, gen- 
erated the interpolated values.) Remember, 
when looking at figure 3, that each interpolation 
arc is a different cubic polynomial, with the 
ensemble of cubics having the unique property 
that they are common and tangent at the base 
points. 
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SCALE OF X 


FIGURE 3. — Four-point inlerpolofion and integration. 

In summary, four-point interpolation estab- 
lishes a form-free ensemble with two orders of 
continuity. Generally, the interpolation locus will 
pass through the base points without undesirable 
oscillation. In some situations, with an abrupt 
change in slope of the base point data, a tendency 
to oscillate may develop. 

Four-point interpolation can be reduced to 
coefficients that operate directly on the base 
points, just as the two-point system was re- 
duced. Consider the special case of 1^1=1 and 
752=753= YB4=0. Then, from equations 33-36 
or equation 37, a=2, 6=-4, c=5/2, and d= 
-1/2, and the interpolation cubic is 

2/o=2-4z-1-{6/2)«2-(1/2)2». (38) 

Equation 38 can be evaluated for selected 
values of z, producing a coefficient that gives the 
unit contribution of the base point 7B1 to the 
interpolated value. Since 7B1 is at a z-value of 
zero, call this coefficient CO. 

'se the special case 7B2=1 and 7B1 
'‘?4=0. The interpolation cubic from 
3-36 is 

yi=8+(19/2)2-722-(-(3/2)2'''. (39) 

Evaluate this for the same selected 2 -values 
used with equation 38. Call this set of coefficients 


Cl, since 752 is at a z-value of 1. 

Similarly, evaluate C2 and C3. An interpolated 
value is the summation of the products of the unit 
operators with any real set of 751 through 754. 
This product summation is 

j/=?/o+vi +^2+vt=coyBi-f Cl- m+c2-yB3 

-t-C3.yB4 (40) 

Table A-1 (appendix A) is a listing of the four 
interpolating coefficients for each one-hundredth 
of the interpolating arc from z—1 to 2=2. These 
coefficients were generated by program A-2. 

The area under any function is given by the 
integral of the function across a range of x (or z). 
The integral of the interpolation cubic of the 
four-point method is given in equation 41. Inte- 
gration, just as interpolation, must be limited to 
the central interval of the four-point set: 

I* j,dz=:a2+(6/2)2H(c/3)23-f(d/4)2‘+C=A. (41) 

The constant of integration, C, can be evalu- 
ated by noting that the area under the curve, A, 
is zero when 2 = 1 ; that is, at the left end of the 
interpolation arc. For this value of 2 , and with 
equations 33-36, C is given by 

C=(l/24)(-17.m+6-K52-19-yB3-f7-yB4). (42) 

The area under any segment of the interpoia- 
tion arc can now be calculated with equations 
33-36, 41, and 42. It is also possible to reduce 
four-point integration to a coefficient-multiplier 
scheme analagous to equation 40 for interpola- 
tions. Table A-2 is a listing of such coefficients, 
also generated by program A-2 (appendix A). 
Now, note that four-point integration is also 
included in program A-1. The area under the arcs 
of the interpolation cubic are plotted in figure 3. 
Remember that each arc starts with zero area. 
The total area to the left of any point is thus the 
sum of the partial area in the arc plus all full arcs 
to the left. 

The coefficients in table A-2 are based on a 
scale variate z with the base points at unit 
intervals of 2 . If the base points are at a scale s in 
the original data, given by s=A:2+Zo, then 
ds=kdz. An integral ydz can be changed to an 
integral yds^^kydz by simple change of variable. 
This change of variable is included in program 
A-1, but must be performed as a separate 
operation following use of table A-2. 
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Six-Point Interpolation 


The development of the method of six-point 
interpolation has been given in detail previously 
(Snyder 1967a). The geometrical and analytical 
concepts are analagous to those used in develop- 
ment of the four-point method in the previous 
section. 

Note briefly that six-point interpolation re- 
quires functions two orders higher than four- 
point interpolation, which in turn was two orders 
higher than the functions of two-point interpola- 
tion. Specifically, six-point interpolation is based 
on a sliding, transhaping, fourth-power polyno- 
mial, two orders up from the sliding parabola. 
Each arc of the interpolation locus is given by a 
fifth-power polynomial. 

Six mathematical restraints are required. In 
addition to the four used for four-point interpola- 
tion, commonality of function and first deriva- 
tive, commonality of the second derivative with 
the second derivative of the linearly transhaping, 
fourth-power polynomial is required at the ends 
of the interpolation arcs. The junctures of the 
arcs at the base points are mathematically 
continuous and have continuous first and second 
derivatives. 

Details of the solution of the coefficients of the 
six-point method are not repeated here. The 
solution for the coefficients in equation 1 is 
conveniently expressed by the matrix-vector 
equation 


(1/24) 


432 -2040 4080 -4080 2040 -408 


rao 


a 

-918 4436 -8762 8712 -4346 868 


YBl 


b 

765 -3754 7414 -7356 3661 -730 


Ym 


c 

-313 1651 -3078 3068 -1621 303 


YB3 


d 

63 -314 626 -624 311 -62 


YBA 


e 

-6 26 -60 60 -25 5 


ras 


f 


(43) 


In the previous section on four-point interpola- 
tion, it was pointed out that the method could be 
extended to integration of experimental data. An 
advantage of the six-point system is that it can be 
converted to differentiation, leading to deter- 
mination of smooth rates of change of experi- 
mental data. The rate of change of values in 
equation 1 is given directly by equation 2. With 
coefficients b through / given by equation 43, 
rates of change can be computed by direct 
substitution in equation 2. An advantage of this 
method of interpolative differentiation is that the 
first derivatives are smoothly continuous be- 
cause the second derivatives are continuous. 


Figure 4 shows six-point interpolation through 
and the calculated first derivative of the same 
data as in figure 3. The interpolative values and 
derivatives were computed with program A-3 
(appendix A). 

Six-point interpolation and differentiation can 
be reduced to multiplying coefficients in a man- 
ner similar to interpolation and integration with 
coefficients as presented in the four-point 
method. Interpolation coefficients for each one- 
hundredth of the interval of the interpolation arc 
from z=2 to z=3 are given in table A-3 (appen- 
dix A). These coefficients were generated by 
program A-4. Table A-4 gives coefficients for use 
in differentiation. Calculation of these coeffi- 
cients is included in program A-4. 

It should be noted that slopes given by the 
differentiation coefficients of table A-4 are in 
terms of the scale variate z. Given that the scale 
of original data, s, is related to z by s~kz+Zo, 
then ds=kdz. A sXopa, dy/dz can be changed to a 
slope of original data by simple change of 
variable, yielding dy/ds~{l/k) dy/dz. This scale 
conversion is included in program A-3, but must 



FIGURE 4, — Six-point interpolation and differentiation. 





be performed as a separate operation following 
use of table A-4. 

SMOOTHING OF 
EXPERIMENTAL DATA 

Smoothing of experimental data is essentially 
the application of the method of least squares to 
effect a transformation from recorded data to a 
set of base points for interpolation. Data may be 
recorded at nonuniform intervals. Transforma- 
tion produces base points at uniform intervals. 
Data may include mensuration errors. These 
errors may be smoothed by averaging during the 
transformation. Successful averaging requires 
that base points be fewer in number than data 
points and also requires a configuration of data 
points conducive to production of good averages. 

The location and spacing of the base points is a 
problem in design. No hard rules can be stated. 
The heart of the problem is the amount of 
averaging that is desired. With a goodly amount 
of data in a compact mass, a high degree of 
averaging can be used. The base points would 
likely be one-fifth to one-tenth of the number of 
data points. With limited data, strung out in 
range, the base points would likely be about 
one-half the number of data points, resulting in a 
limited degree of smoothing. 

In the following examples, a high degree of 
smoothing will be evident. The placement of 
points is usually dictated by the particular 
problem. A base point would usually be placed at 
an interface or boundary, such as the surface of 
ground or water. In time series one point would 
usually be placed at time zero. Spacing of the 
remaining base points through the data range is 
obviously dictated by the degree of averaging 
and hence the number of points. 

Two-Poinf Smoothing 

Two-point smoothing is identical to evaluation 
of connected linear segments (Snyder 1967b). 
Equation 14 can be generalized as 

y=:C0-YB0+Gl-YBl. (44) 

In equation 44, the meaning is limited to one 
interpolation interval. If several adjacent inter- 
vals are considered, equations 45-47 result: 

i/,=co-yBO-f-ci']rai (45) 

y2= CO-rai-hCl-m (46) 

yj= co-ym-t-ci-ysa (47) 


In equation 45, y is any value lying on the line 
from YBO to YBl. In equation 46, y., is any value 
lying on the line from YBl to YB2, and so on. The 
interpolation coefficients CO and Cl are known 
functions of scale variate z, as in table 1. 

Assume that experimental data consist of 
simultaneous coordinate values of x and y for a 
series of observations. Then, consider a series of 
base points to be substituted for the data points. 
It is necessary to specify the location of the base 
points within the x-scale of the data. Normally, 
one has some idea of the number of base points 
desired to span the data in the x-scale. Setting 
the location of the base points in x-scale automat- 
ically sets a z-scale value for each data point since 
the base points are always spaced at unit interval 
in z. 

Referring now to equations 45-47, note that 
they are linear in form. Note also that given 
experimental values of y and given a computed 
z-variate for each observation, there is an equiv- 
alent of a series of observational equations as 
used in multiple regression. CO and Cl take the 
position of x„ and x, in regression, and YBO and 
YBl take the position of regression coefficients, 
conventionally designed 6o and bj. There is no 
regression intercept normally designated a. It is 
a simple matter, with conventional least-squai’es 
methods or programs, to compute YBO and YBl 
as regression coefficients, A transformation from 
experimental data in x-y coordinates to base 
point data, in z-YBl coordinates, is thus accom- 
plished. Note especially that in the configuration 
of equations 45-47, YBl and YB2 must take 
values on two adjacent intervals. The continuity 
of the smoothing operation is thus satisfied with 
no gaps at the base points. 

Table B-1 (appendix B) shows the details for 
organization of 30 data points into format for 
regression transformation, The data points and 
the result of least-squares transformation in 
two-point smoothing are plotted in figure 5. 
While the location of base points at x=0.4, 0.6, 
0,8, and 1.0 is arbitrary, the values of these base 
points are wholly objective and the result of 
least-squares fitting. Note specifically in figure 6 
the overall fit of the linear segments to the rough 
crescent shape of the experimental data. Since 
the two internal base points are each required to 
be common to two segments, the fit of a single 
segment to the data points in that segment is 
secondary. 
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FIGURE 5. — Smoothing by two-point method. 

Four-Point Smoothing 

Four-point smoothing is a simple extension of 
two-point smoothing. Equation 40, modified for 
adjoining arcs, is extended in equations 48-51, 
analagous to the two-point equations, 45-47. 

V^^CQ'YBQ+ChYBl+Ci'Ym-VC^-YBi (48) 

C0-YB2+ChYB2+C2’YBa+C3-YB4 (49) 

2 / 3 = CQ-YB2+ChYB3 

-hC2>yB4-fC3.yB6 (60) 

2/1= co.ys3-fci-ra4 

+C2-ra6-hC3.KB6 (61) 


Values for the coefficients are given in table 
A-1. Remember, in equations 46-47 each base 
point is linked to two adjacent intervals. FSl and 
yB2 in those equations are each linked by CO and 
Cl. Note in equations 48-51 that an internal base 
point such as YB3 is linked to four interpolation 
arcs. The least-squares value of YB3 is depend- 
ent on experimental data in four intervals. 
Exterior base points are dependent on 1, 2, or 3 
intervals. 

Table B-2 (appendix B) shows organization of 
experimental data in format for four-point 
smoothing. The same data are used as for the 
two-point smoothing example, except the a;-scale 
is shifted to the left to put the first base point at 
x=0.4. The 2 :-scale is not affected. Program B-1 
(appendix B) is used to perform the least-squares 
transformation. This program also generates the 
locus of interpolated values which runs through 
the derived base points and generates the inter- 
grals of the interpolation arcs. 

Figure 6 shows the results of the four-point 
smoothing. An interpolation locus of smoothly 
joined arcs is generated that runs continuously 
through the data. Comparison of figures 6 and 6 



FIGURE 6. — Smoothing by four-point method. 


shows that the two-point and four-point methods 
produce base points not greatly different. The 
main difference is the creation of arcs by the 
four-point method, smoothly joined with angles 
eliminated. Values of the derived base points can 
be compared in table 2. 

The integral of the interpolation arcs is also 
plotted in figure 6. This integral can be thought 
of as the smoothed summation of the expei’i- 
mental data. 

Six-Point Smoothing 

Six-point smoothing is a simple extension of 
methodology already presented. Equations 52-65 
are analogous to equations 45-47 and 48-61. 

2 /, =C0-yB0-|-Cl-yBl-t-C2-KB2 

-fC3-ra3-fC4.KS4-fC6'yB6 (62) 

2/ =C0- KBIH-Cl- m-fC2' m 

+C3>YBA+CA-YB5+C5-YBe (63) 

y.=C0-YB2+ChyB3+C2-YB4 

+C3-yB^+CA-YB6+C5‘yB7 (64) 

2/j=C0-ra3-(-Cl'yB4-t-C2.YB5 

-f-C3-ra6t-C4.yB7-fC6'yB8 (65) 

In equations 52-55, note that an interior base 
point will depend upon values of experimental 


Table 2. — Comparison of derived base points 


Two -Point 

Foiir-Pomt 


Sjx.Pomt 

X 

YBi 

X 

YBi 

X 

YBi 





0.4 

-4.14772 



0.4 

-1.72742 

.6 

-1,65936 

0.4 

1.62418 

.6 

1.42889 

.8 

1.417602 

.6 

3.76141 

.8 

3.69534 

1,0 

3.69252 


3.81651 

1.0 

3.67674 

1.2 

3.68015 

1.2 

2.67954 

1.2 

2.44932 

1.4 

2.46502 



1.4 

- .42236 

1.6 

-2.29737 

♦ ‘ 


.... 



1.8 

-20.17608 
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data in six adjacent intervals by linkage through 
the coefficients. Exterior base points may de- 
pend on five or fewer intervals of data. 

Table B-3 (appendix B) shows organization of 
experimental data in format for six-point smooth- 
ing. Program B-2 is used to perform the least- 
squares transformation. Following transforma- 
tion of the experimental data to a smoothed set of 
base points, program B-2 computes a locus of 
interpolated values and a locus of differentials of 
the interpolation arcs. Figure 7 is a plot of the 
experimental data and derived interpolation and 
differentiation loci. 



Figure 7 contains the same experimental data 
as figures 5 and 6. It can be seen that the 
interpolation arcs of the four-point and six-point 
methods are almost identical. Numerical com- 
parison of the derived base values is shown in 
table 2. 

Figure 7 also shows the differentiation of the 
interpolation arcs. This locus of values can be 
regarded as the smoothed rate of change of the 
experimental data. Note the smoothness of this 
locus through the base points and the near 
constant value of slope from about *=1.14 to 
«= 1 . 2 . 

ATTRIBUTES OF 
-SLIDING POLYNOMIALS 

'’.and 7 demonstrate that various 
lethods can be used to smooth 
ata. Smoothed values, whether 
>-point, four-point, or six-point 


interpolation loci, are not extremely different. 
Other methods would also produce acceptable 
interpolation loci. Kimball’s (1974) method would 
require a means to separate the spectral portions 
of noise and information in the data. 

The identifiable advantages of sliding polyno- 
mials are in the secondary characteristics of their 
loci. The four-point and six-point methods are 
smooth because their derivatives are continuous 
at the juncture points of the interpolation arcs. 
The four-point method, with a continuous first 
derivative, produces an extremely smooth locus 
of integral values, since integration itself is a 
smoothing process. The six-point method, with 
continuous first and second derivatives, can gen- 
erate a smooth rate of change in experimental 
data. DuChateau et al. (1972) also used six points 
in evaluation. In their method, however, para- 
bolic splines are required to be continuous only in 
the first derivative. Therefore, these splines 
cannot be expected to develop the smoothly 
continuous rates of change of experimental data. 
Such continuous smoothness is dependent upon a 
continuous second derivative. 

Sliding polynomials and their integrals and 
derivatives can be derived from experimental 
data by conventional least-squares techniques. 
The interpolation coefficients provide exact error 
weights for each experimental data point, re- 
gardless of the interpolation interval in which the 
data point is located. The structure of equations 
52-55, for example, shows that the evaluation 
procedure is similar to the use of smoothing 
convolutes as in Savitzky and Golay (1964), 
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APPENDIX A.— COEFFICIENTS AND COMPUTER PROGRAMS FOR 

INTERPOLATION 

table A'l . — Coefficients for four-point interpolation 


1.00 

0.00 000000 

1.00000000 

. 00 06000 0 

0.00000000 

wsaam 

-.00490050 

-.00960400 

.99975450 

.99901200 

. 00519850 
.01078800 

00004950 

-.00019600 

1.03 

1.04 

-.01411350 

-.01843200 

.99779050 
.9960960 0 . 

.01675950 

.02310400 

-.00043650 

-.00076800 


1.05 

1.06 

-.02256250 

-.02650800 

.99393750 

.99132400 

.02981250 

.03667600 

-.00118750 

-.00169200 

m 

-.03027150 

-.03385600 

. 98826450 
.98476800 

.0 442'8550 
.05203206 

-.00227850 

-.00294400 

1.09 

-.03726450 

.98084350 

.06010650 

-.00368550 

1.10 

-.04050000 

.97650000 

.06650000 

-.00450000 

1.11 

1.12 

-.04356550 

-.04646400 

.97174650 

.96659200 

.07720350 

.06620800 

-.00538450 

-.00633600 

1.13 

1*14 

-.04919850 

-.05177200 

.96104550 

.95511600 

.09550450 

.10508400 

-.00735150 

-.00842800 


1.15 

1.16 

-.05418750 

-.05644800 

.94881250 

.94214400 

.11493750 

.12505600 

-.00956250 

-.01075200 

1.17 

1.16 

-.05655650 

-.06051600 

.93511950 

.92774800 

.13543050 

.14605200 

-.01199350 

-.01328400 

1. 19 

-.06232950 

.92003850 

.15691150 

-.01462050 

1.20 

-.06400000 

.91200000 

.16800000 

-.01600000 

1.21 

1.22 

-.06553050 

-.06692400 

.90364150 

.89497200 

.17930850 

.19062800 

-. 017419'50 
-.01887600 

1.23 

1.24 

-.06818350 

-.06931200 

.88600050 

.87673600 

.20254950 

.21446400 

-.02036650 

-.02188600 


1.25 

1.26 

-.07031250 

-.07116800 

.86718750 

.65736400 

.22656250 

.23883600 

-.02343750 

-.02501200 

1.27 

1.26 

-.07194150 

-.07257600 

.84727450 

.83692800 

.25127550 

.26387200 

-. 02660850 
-.02822400 

1.29 

-.07309450 

.62633350 

.27661650 

-♦02985550 

1.30 

-.07350000 

.81550060 

.28950000 

-.03150000 

1*31 

1*32 

-.07379550 

-.07398400 

•80443650 

.79315200 

*30 251350 
.31564800 

-.03315450 
-. 0348160.0 

1.33 

1.34 

-.07406850 

-.07405200 

.78165550 

.76995600 

.32889450 

.34224400 

-.03648150 
-. 03814800 







Table k-l.— Coefficients for four-point interpolation —Continued 


1.35 

1.36 

-.07393750 

-.07372800 

.75806250 

.74598400 

.35568750 

.36921600 

-.03981250 

-.04147200 

1.37 

1.38 

-.07342650 

-.07303600 

.7337295 0 
.72130800 

.38282050 

.39649200 

-.04312350 

-.04476400 

1.39 

-.07255950 

.70672850 

.41022150 

-.04639050 

1.40 

-.07200000 

.69600000 

,42400000 

-.04800000 

1.41 

1.42 

-.07136050 

-.07064400 

.68313150 

.67013200 

, 43781850 " 

. 45166800 

-.04958950 

-.05115600 

1.43 

1.44 

-.06985350 

-.06899200 

.65701050 

.64377600 

,4655 395 0 " 

, 479*^2400 

-.05269650 

-.05420800 

1.45 

1.46 

1.47 

1.48 

-.06806250 

-.06706800 

-.06601150 

-.06489600 

.63043750 
. 61700400 
.60348450 
.58988800 

.49331250 
.50 719600 
.52106550 
.53491200 

-.05568750 
-.05713200 
-.05853850 
-. 05990400 

1.49 

-.06372450 

,57622350 

,54872650 

-.06122550 

1.60 

-.06250300 

.56250000 

.56250000 

-.06250000 

1.51 

1.52 

-.06122550 

-.05990400 

.54872650 
. 53491200 

.57622350 

.58988800 

-.06372450 

-.06489600 

1.53 

1.54 

-.05853850 

-.05713200 

,5 210655 0 
,50719600 

.60348450 

•61700400 

06601150 

-.06706800 


1.55 

1.56 

-.05568750 

-.05420800 

.49331250 

.47942400 

.63043750 

.64377600 

-. 06806250 
-.06699200 

1.57 

1.58 

-.05269650 

-.05115600 

.4655395 0 
.45166800 

.65701050 

.67013200 

-.06985350 
- • 07064400 

1.59 

-.04958950 

,43781850 

.68313150 

-.07136050 

1.60 

-.04800000 

. 42 hOOOOO 

.69600000 

-.07200000 

1.61 

1.62 

-.04639050 

-.04476400 

.41022150 

,39649200 

.70872850 

.72130800 

-. 07 255950 
-.07303600 

1.63 

1.64 

-.04312350 

-.04147200 

.38282050 

.36921600 

.73372950 

.74598400 

-. 07342650 
-. 07372800 


1.65 

1.66 

-.03981250 

-.03814800 

.35566750 
* 34224400 

.75806250 

.76995600 

-. 07 393750 
-.07405200 

1.67 

1.68 

-.03648150 

-.03481600 

,32889450 
. 31564800 

.78165550 

.79315200 

-. 07406850 
-.07398400 

1.69 

-.03315450 

.30251350 

.80443650 

-.07 379550 



Table A-l . — Coefficients for four-point inteTpolation--Coni\r\\xe^ 


1.70 

>.03150000 

. 28950000 

. 81550000 

-.07350000 

1.71 

1.72 

-.02985550 

-.02822400 

.27661650 
. 26387200 

.82633350 
. 83692600 

-.07309450 

-.07257600 

1.73 

1. 74 

-.02660950 

-.02501200 

.25127550 
. 23883600 

.84727450 
. 85736400 

-.07194150 

-.07118800 


1.75 

1.76 

-.02343750 

-.02188800 

.22656250 
, 21446400 

.86716750 
. 87673600 

-.07031250 

-.06931200 

1.77 

1,76 

- ,02036650 
- ,01887600 

.20254950 
. 19082800 

♦86600050 
, 89497200 

-.06818350 

-.06692400 

1.79 

-.01741950 

.17930850 

.90364150 

-.06553050 

1.60 

-.01600000 

.16600000 

.912G0000 

-.06400000 

1.81 

1.82 

-. 01462050 
-.01328400 

,15691150 
. 14605200 

.92003850 
. 92774800 

-.06232950 

-.06051600 

1.63 

1.64 

-.01199350 
- ,01075200 

. 13543050 
, 12505600 

. 93511950 
. 9421440 0 

-.05855650“ 

-.05644800 


1.86 

1.86 

-.00956250 

-.00842800 

.11493750 

.10508400 

.94881250 

.95611600 

-.05418750 

-.05177200 

1.87 

1.86 

-.00735150 

-.00633600 

.09550460 

.08620800 

.96104650 
, 96659200 

-.04919850“ 

04646400 

1.89 

-.0 05 38 450“ 

,07720350 

.97174650 

-.04356550“ 

1.90 

-.00450000 

.06850000 

.97650000 

-.04050000 

1.91 

1.92 

-.00368550 

-.00294400 

.06010650 
. 05203200 

.98084350 

,98476800 

-.03726450 

-.03385600 

1.93 

1,94' 

- ,00 227850 
-.001 69 200 

♦04428550 

•03687600 

,98 82 645 0 
.99132400 

-.03027150 
. -,02650800 

1.95 

1.96 

1.97 

1.98 

-.00118750 
- ,00076800 
-.00043650 
-.00019600 

.02981250 

.02310400 

.01675950 

•01078800 

.99393750 

.99609600 

.99779050 

.99901200 

02256250 

-.01843200 

01411350 

-.00960400 

1.99 

-.00004950 

.00519860 

.99975150 

-.00490050 

2.00 

0,00000000 

-.00000000 

1. 00 OUOOOO 

0. OQOOOOOO 



Table A'2.— Coefficients for four-point integration 


X 


cr 

C2 

C3 



} 


1 . ‘ 

1*00 

>.00000000 

-* OOQOOOliO 

O.OOOGQOOO 

.00000000 

1.01 

1.02 

-.00002467 

-.00009735 

.00999917 

.01999339 

.60002566 
, 00010527 

-.60000017 

-.00000131 

1.03 

1.04 

-.00021610 

-.00037899 

.02997780 

.03994763 

. 00024270 
,00044171 

-. 00000440 
-.00001035 


1.05 

1.06 

-.00058411 

-.00082962 

.04989818 

.05982486 

.00070599 

.00103914 

-. 00002005 
-.00003438 

1.07 

1.06 

-.00111367 

-.00143445 

.06972317 

.07958869 

.00 144466 
.00192597 

-.00005417 

-.00008021 

1.09 

-.00179020 

.08941710 

*00248640 

-, 00011330 

1.10 , 

-.00217917 

,09920417 

.00312917 

-.00015417 

1.11 

1.12 

-.00259963 

-.00304992 

♦10894574 

.11863776 

.00385743 

,00467424 

-.00020353 

-.00026208 

1.13 

1.14 

-.00352837 

-.00403335 

.12827627 

,13785739 

.00558256 

.00658527 

-,00033047 

-.00040931 


1.15 

1.16 

-.00456328 

-.00511659 

.14737734 

,15683243 

,00768516 

.00888491 

-.00049922 

-.00060075 

1.17 

1.18 

-.00569173 

-.00628722 

.16621904 

.17553366 

.01018713 

.01159434 

-.00071443 

-.00064078 

1.19 

-.00690157 

.18477287 

.01310896 

-.00098027 

1.20 

-.00753333 

,19393333 

.01473333 

-.00113333 

1*21 

1.22 

-.00818110 

-.00884349 

,20301180 

.21200513 

.01646970 

.01832021 

-.00130040 

-.00148185 

1.23 

1.24 

-.00'951913 

-.01020672 

.22091024 

,22972416 

.02028693 

.02237184 

00167803 

-.00188928 


1.25 

1.26 

-.01090495 

-.01161255 

, 23844401 
,24706699 

.02457682 

.02690367 

-. 00211589 
-.00235811 

1.27 

1.28 

-.01232830 

-.D'1305099 

,25559040 

,26H01i63 

.02935410 

,03192971 

-. 00 261620 
-.00289035 

1.29 

-.01377943 

,27232814 

.03463203 

-. 00318073 

1*30 

-.01451250 

.28053750 

.03746250 

-.00348750 

1.31 

1.32 

-.01524907 

-.01598805 

,28863737 

,29662549 

. 04042246 
,04351317 

-. 00381077 
-.00415061 

1.33 

1.34 

-.01672840 

-.01746909 

,30449970 

,31225793 

,04673580 

.05009141 

-.00450710 
-. 00488025 
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Table A-2.— Coefficients for four-point integration — Continued 


1.35 

1.36 

-.01820911 

-.01894752 

.31989818 

.32741856 

.05358099 

.05720544 

-.00527005 

-.00567648 

1.37 
1. 38 

-.01968337 

-.02041576 

.33481727 

.34209259 

.06096556 

.06486207 

-.00609947 

-.00653891 

1.39 

-.02114380 

.34924290 

.06889560 

-.00699470 

1.40 

-.02186667 

. 35626667 

.07306667 

-.00746667 

1.41 

1.42 

-.02258353 

-.02329362 

.36316244 
.36992 886 

.07737573 

.08182314 

-.00795463 

-.00845830 

1.43 
1. 44 

-.02399617 

-.02469045 

.37656467 

.38306869 

.06640916 
. 09113397 

-.00 897767 
-.00951221 

1.45 

1.46 

1.47 

1.48 

-.02537578 

-.02605149 

-.02671693 

-.02737152 

.38943984 

.39567713 

.40177964 
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-.01120423 
-.01179648 
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-.02801467 
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-.01302083 

1.51 
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-.0 3'320 829 
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Table A-2. — Coefficients for four-point integration — Continued 


\.70 

-.03817917 

.50420417 

.26112917 

-.02715417 

1.71 

1.72 

-.03848593 

-.03877632 

.50703464 

.50973696 

.26933853 

.27765504 

-.02768723 

-.02861568 

1.73 

1.74 
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1.75 

1.76 
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“•.03977739 
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-.03076172 

-.03145995 

1.77 

1.78 

-.03998863 
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. 52334646 

.32075643 

.32966154 

-.03214753 

03262318 

1 . 7 <} 
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.33865486 

-.03348557 

1.60 

-.04053333 

.52693333 

.34773333 

-. 03413333 
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1.82 

-.04068640 
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-.03476510 

-.03537945 

1.63 

1*84 
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-.04106592 
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.37544763 

.38483424 

-.03597493 

-.03655008 


1.85 

1.66 

-.04116745 

-.04125735 

.53398151 

.53508139 

.39428932 

.40380927 

-.03710339 
-. 03763331 

1.87 

1.86 

-.04133620 
- • 04140 459 

.53608410 

.53699243 

.41339040 
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-.03813830 
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1.89 

-.04146313 
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1.90 
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-.04083705 
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1.96 

1.97 
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.54156139 

.49176849 
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.51168886 

.52167327 

-. 04108255 
-.04128768 
-.04145057 
-.04156931 

i . 99 

-.04166650 

.5416410 0 

.53166750 

-.04164200 

2.00 

-.04166667 

.54166667 

.54166667 

-. 04166667 
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Table A-3. Coefficients for six-point interpolation — Continued 
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Table A-3. — Coefficients for six-pmnt interpolation — Continued 
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Table A-4. — Coefficients for six-point differentiation 
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A -4. — Coefficients for six-point differentiation — Continued 
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Table A-4. — Coefficients for six-point differentiation — Continued 
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Table A-4. — Coefficients for six-point differentiation — Continued 
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EXAMPLES USING TABLES A-1 AND A-2 


EXAMPLES USING TABLES A-3 AND A-4 


Example 1, — Interpolate a value of Y at x=0.07 
Given the table of base points 


X 

BYt 

-0.05 

1.00 

.00 

1.15 

.05 

1.25 

.10 

1.25 

.15 

1.10 

.20 

.75 

.25 

.15 


2=1 + (0.07-0,05)|10.10-0.05)=1.4 

From table A-1: 


(C;0= -0.0721x1.15= 

-0.0828 

(Cl= .6961x1.25= 

.8700 

((72= .4241x1.25= 

.5300 

(<73= -.0481x1.10= 

-.0518 

y 0B^ 

1.2644 


Example 2. — Compute the integral from x=0.00 
to x = 0.07 


From 0. 00 to 0. 05 
From table A-2: 


(C0= 

-0.041671x1.00= 

-0.0417 

(Cl= 

.541671x1.15= 

.6229 

(C2= 

.541671x1.25= 

.6771 

(C3= 

-.041671x1.25= 

-.0521 



1.2062 

From 0.05 to 0.07 

From table A-2: 


(C0= 

-0.021871x1.15= 

-0.0252 

(Cl= 

.356271x1.25= 

.4453 

(C2= 

.073071x1.25= 

.0913 

(C3= 

-.007471x1.10= 

-.0082 



.5032 

/'0 07 

-^.00 =11.20624-0.50321x0.05=0.08547 


Example 3. — Given the same base points as in 
example 1 interpolate at x=0.07. 

From table A-3: 

(C0=0.0144)xl.00= 0.0144 
(C'l=-.1136)xl.l5= -.1306 
(C2= .72641x1.25= .9080 

(C3= .43841x1.26= .5480 

(C4=-.07361xl.l0= -.0810 
{C5= .00801 X .75= .0060 

(107 = 1.2648 


Example 4. — Interpolate for the differential of y 
at x=0.07 

From table A-4: 

ini=0.018001xl.00= 0.0180 

(01= .097331x1.15= .1119 

I r2= 1.322671x1.25= 1.6533 

(GO =1.464001x1.25= 1.8300 

(C4= .259331x1.10= -.2853 
(C5= .03867} X .75= .0290 

1)143 


=(O.OM3/0.05}=0.286 


NOTES FOR PROGRAMS A-1, A-2, A-3, AND A-4 

TAIKIl Input of tabular values, or base 
points, on which interpolation is 
based. 

YINTIH The interpolated values. 

YAlNTdl Interpolated integral. 

YPINTdl Interpolated derivative. 

INTFORl 1 Subroutine for four-point 
interpolation. 

INTSIXl 1 Subroutine for six-point 
interpolation. 

dd) Coefficients of equations 37 and 43. 

IGOEFd) The matrices of equations 37 and 
43. 

YAOdl Coefficients of equation 42. 

Xldl The index z of tables A-1, A-2, A-3, 

and A-4. 
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PROGRAM A-1. — Four-point interpolation and integration 


DIMENSION YINT(200),YAINT(200) 
COMMON/XX/TAP(100) 

READ(5,1000) {TA3(I) ,1=3,22) 

1000 FORMAT{20F4,2) 

TAP (1) = 0. 1 
TAP(2)=0. 1 
X=0.3 

MS=-1 

no 1001 1=1,151 
X=X+0.01 

CALL INTFOR(X,Y,YA,MS) 

YINT(I)=Y 
YAINT(I) = YA 

1001 CONTINUE 

WRITE (6, 100?) { I, YI NT (I) ,1=1,151) 

1002 F0RMAT(» INTERPOLATED FUNCTION^/ { 8 { 15, FIO . 4) ) ) 
WRITE (6, 100 3) (I,YAINT(I) ,1=1,151) 

10 0 3 FORMAK^ INTERPOLATED INTEGRAL^/ (8 (15, F 10 . 4 ) ) ) 
STOP 
END 


SUPROUTINE INTF0R(X,Y,YA,MS) 
dimension C(4) ,IC0EF (4,4) 

COMMON/XX/TAP(100 ) 

DATA IC0EF/4,-R,5,-l,"6,19,-14,3,6,-16, 13,“3,-2,5,-4,l/ 
AX=((X-TAB(1))/TAB(2)) 

M=INT(AX) 

FX=AX-FLOAT(M)+1.0 
IF(M.EQ.MS) go TQ J02 
DO 100 1=1 ,4 
C{I)=0.0 
00 100 J=l,4 

10 0 C (I)=C(I) +ICOEF(I, J)^TA8(M+J+1) 

CA=C(l)+O(2)/2.a + C(3)/3.0t-C(4) /4.0 
MS = M 

102 Y = C(4) 

Y A=C(4) /4, 0 
DO 101 K=l,3 
Y=Y»FX+C(4-K) 

101 YA=YA»FXt-C(4-K)/(4-K) 

Y=Y/2,0 

YA=(YA»FX-CA)*TAR(2)/2,0 

RETURN 

END 
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PROGRAM A-2. — Generation of four-point coefficients 


iiMENSiON YiNtdQl ,4) ,'Y4INT (101,4) ,IC0EF(4,4) ,XI (101) , YA0(4) 
J f" / 4," 8 , 5 , - 1 6, 19_,-1 4 , 3 ,6 , - 16 , 13 3 , -2 , 5 , -4 , 1 / 

YAO(l) =-17.0/24.0 
J A0(^)_=5. C/24^0 

Y AO ( 3 ) = - 19. 07 2 4 i'O 
YAO (4) =7. 0/24. 0 

"oo' iOGO J = i',i6i ' " 

X=FLOAr (J)/ 10G.J »C.99 

xi(j)=x 

00 10^1 ^J.,4 _ 

YiNY(j,i) =ic6ef(4', I)' ■ 

00 1002 K = l_,3 

1 0'O 2 Y INt ('jTiT'= Y fNT"( jriT*X + IC0EF (4- K ,1) 

YINT( J,I) =YINT <J, I )/2.0 

Y A I NT ( J , I ) = I COE F ( 4 ,T ) >4 . 0 
□ O_lJ!0 J J< = 1,3 

10 0 3 ■'YAINf'f I ) = { Y AINf ( J, I ) *X ) + { FL OAT ( ICOEF ( 4- K , I ) ) / (4- K) ) 
YAINT(J,I) = (YAINT ( J , I ) /2 ♦ 0 ) »X-t-YA0 ( I ) 

lOOi “OONTINUE" ' 

1000 CONTINUE 

do 10 0 5 J = l,'i01 

IF(J.E3.1 .OR.J.EI^. 51) GO TO 1006 
GO TO 1007 

1006 WRITE (6_, ID 082 
lb0 8"FdRMAT(>r») 

1007 IF(MOD (^,_52.. EGlil2_G^ TO ^ 

GO TO lOib 

100 9 WRITEiSilCUL 

ioil FORMAT (♦0+) 

10 10 WRI TE(6,1C12 ) (XKJ), < YINT(J,I) ,1=1,4)) 
i Oia'KORMATf* ♦,(35)<,F5i2“,"4Fi2.8y) 

10 05 conjxnue _ 

DO 1020 J=l,101 

1 f" < . L 1 .OR . J.EO . 51) G O TO _1 0 13 
■ 'Gb”TO iT14 

1013 WRITE(6,10C8) 

1 0 1 4"l F ( MOD ( J, 5TrEO . £r“GO TO" IOTG ' ' 

GO TO 101 7 _ 

lO'ie WRltF(6,l'0lld 

1017 NRITE(6,1Q12) (XKJ) , (Y/U_NT(J,I) ,1=1,4) ) 

102 0' CONTINUE 

STOP 
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PROGRAM A-3.- — Six-point interpolation and differentiation 


DIMENSION YINT{200)>YPINT(200) 

COMMON/XX/TAPttOO) 

READ(5,lQOOt <TA0(I),I=3t22) 
lOOD P0RMAT<2QF4.2) 

TA3(1>=0. 1 

TA0(2)=O.l 

X = 0,3 

MS=-1 

DO 1001 1=1,151 
X=XfO. 01 

CALL INTSIX(X,Y,YP,MS) 

YIMT(I)=Y 

YPINT(I) = YP ■ ■■■ 

1001 CONTINUE 

HRITE(6,10 0 2) (ItYINTd) ,I=1,1R1) 

1002 FORMATS INTERPOLATED FUNCTION^/ (ft < 15, F 10 . 4) ) ) 
WRITE(5,i003) <I,YPINT(I) ,1=1,151) 

1003 FORNAT(» INTERPOLATED DE RIVATI VE^/ C 8 (15 , FiO • 4) 1 ) 
STOP 

END 


SUBROUTINE INTSIX (X,Y ,YP ,MS) 

DIMENSION 0(6) ,IC0EF(6,6) 

COMMON/XX/TAPdDO ) 

DATA ICOEF/432, -Pi 8, 765, -313, 6 3, -5, “20 4 0,44 36, -375 4, 1551, -31 4, 25, 
14 080, -875 2,7414, -3 078,62 6, -50, -408 0,871 2, -7356, 3 05 8, -624, 50, 2 040, 
2-4346,3661 ,-1521, 311,-25,-408,868,-730, 303,-62,5/ 

AX={(y-TAB<l) )/TAe(2) ) 

M=INT( AX1 

FX=AX-FL0AT(M)+2. 0 
IF(M.EQ.MSI GO TO 102 
DO 100 1=1,6 
C{I)=0.0 
DO 100 J=l,6 

too C {I)=C(T) + ICOEF (I, J)*TAB{M+J) 

HS = M 

102 Y=C{6) »FX+C(5) 

YP=5.0*C(6) 

00 101 K=l,4 
Y=Y^FX+C(5-K> 

101 YP=YD‘»FX+<5-K)*C<6-K> 

Y=Y/24, 0 

YP=YP/24. O/TA0{2) 

RETURN 

END 
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PROGRAM A-4. — Generation of six-point coefficients 


OIMENSXOM YINT(101,6) , YPINT (IQI ,6) ,IC0EF(6>6> ♦XT (101) 
£lATA_I£;OHF/ii32>-9ia4 765,-313,6 3,-5,-20 40, 4436 |-3754,1551,-314,?5-, 
14 0 8C ,-«75 2,7414,-3 07 8,626,-50 ,-40 8P , 871 2 , -735b, 3 05 8,-624, 50 , 20 40 , 
2-43-46, 366.±,- 1521,311,-25, -40 8, 868,-730, 30 3,-6 2, 5/ 

00 1003 J=l,101 
_X=F.Li3AT (J 1 ymo^ Ati ,.99 
XI(J)=X 

- D.O .IQAl X=l,6 

YINT(J,I)=TC0EF(6,I) 

- OCL 1002 K=l,5 . 

LOO 2 YINT(J,I) =YINT(J,I)'^X + TC0EF(6-X,I) 

YIMT (J, I) =YIMT( 3,11/24.0 
YPINT( J,I ) =IG0EF(6 ,1) ^5, 0 
90 1003 K=l,4 

LOO 3 YPINT(J,I>=yPINT(J,I)»X + TCOEF(6-K,I)*(5-K) 
YPIMTLJ,I)=YPINT(J,I)/24.0 
LOOl GONTINU'^ 

100 0 CONTINUE . .. 

no 1005 j=i,iai 

I-F(a-£Q.l.nR. J.sa. 51) GO TO 1006 
CO TO 100 7 
LQ06 WRITE(6,1008) 

1008 FORMflTt'^l*) 

L007 IF.JMOn ( J,5) .EO.l) GO TO 1009 
CO TO lOlC 

1009 WRlTE(6,10il) 

LOll FCRMATC^O'f) 

1010 NRITE (6,1012) (XL(J) , < YINT ( J , I ) , 1=1 ,6) ) 

101? FORMATt^ (35X,F5.2,6F12,8)) 

L 00 5 CONTINUE , ... 

DO 10 20 J=l,101 

IF( J.ET.l .OR,J.E0.51) GO TO 1013 
CO TO 1014 

1013 WRITE(6,1CP8) 

1014 IF(MOO ( J,5) .FO. 1) CO TO 1016 
CO. TO 1017 

1016 WRITE (5,1011) 

IQIT WRITE(6,1012) (XI(J), (YPINT(J,I),I=1,6)) 

1020 CONTINUE 
CTOP 
c’Nn 


ro 



APPENDIX B. — DATA ARRANGEMENT AND COMPUTER PROGRAMS FOR 

SMOOTHING 


Table B -1 . — Data arrangement for two-point smoothing 


Measured Base-point Scale Values of coefficients^ 

data points positions vanate for base points 


X 

y 

in X 

z 

m 

YB2 

YBS 

YBA 

y 



^ 0,4 







0.412 

1.50 


0.06 

0.94 

0.06 



1.50 

.416 

1,69 


.08 

.92 

.08 



1.69 

.458 

2.02 


.29 

,71 

.29 



2.02 

.450 

2.20 


,25 

.75 

.25 



2.20 

.464 

2.50 


.32 

.68 

.32 



2.50 

,500 

2.50 


.50 

• .50 

.50 



2.50 

,500 

2.70 


.50 

.50 

.50 



2.70 

.514 

2,89 


.57 

.42 

,57 



2.89 

.566 

3.25 


.83 

,17 

.83 



3.25 

,574 

3.58 


,87 

.13 

.87 



3.58 

.600 

3.62 

.6 

1.00 


1.00 



3.62 

.624 

3.66 


.12 


.88 

0.12 


3.66 

.626 

3,85 


.13 


.87 

.13 


3.85 

,660 

3.55 


.30 


.70 

.30 


3.55 

,660 

4.04 


,30 


,70 

.30 


4.04 

.686 

3.88 


.43 


,57 

.43 


3.88 

.698 

3.60 


.49 

, 

.51 

.49 


3.60 

.716 

3,98 


.58 


.42 

.58 


3.98 

,750 

3.70 


.75 


.25 

.75 


3.70 

,774 

3.87 

.8 

,87 


.13 

.87 


3,87 

.806 

3.58 


.03 



,97 

.03 

3.58 

.844 

3.69 


.22 



,78 

.22 

3.69 

,860 

3.29 


30 



.70 

.30 

3.29 

.866 

3,56 


,33 



.67 

.33 

3.56 

.886 

3.27 


,43 



.57 

,43 

3,27 

.910 

3.10 


.55 



.45 

.55 

3,10 

.914 

3.23 


.57 



.43 

,57 

3,23 

.956 

2.74 


.78 



.22 

.78 

2.74 

.960 

2.98 


.80 



.20 

.80 

2.98 

.986 

2,55 

1.0 

.93 


• 

.07 

.93 

2.55 





^1 

^2 ' 

^3 

^4 

r 


Conventional regression symbols 


^ Equations 12 and 13. 
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Table B-2. --Data arrangement for four-point smoothing 



1.06 

-0.0265 

0.9913 

0.0369 

-0.0017 

, 

. . 

1.60 

1.08 

-.0339 

.9848 

.0520 

-.0029 


. 

1.69 

1.29 

-.0731 

.8263 

.2766 

-.0299 

, . 


2.02 

1,25 

-.0703 

.8672 

.2266 

-.0234 

, 


2.20 

1.32 

-.0740 

.7932 

,3156 

-.0348 

. . . . 

... 

2.60 

1.60 

-.0626 

.6626 

.5626 

-.0625 

. . 

. 

2.50 

1,50 

-.0625 

.5625 

.5625 

-.0625 


.... 

2.70 

1.57 

-.0527 

.4655 

.6570 

-.0699 



2.89 

1,83 

-.0120 

.1354 

.9351 

-.0586 

, . 


3.25 

1.87 

-.0074 

,0965 

.9610 

-.0492 




3.58 

2.00 

.0000 

.0000 

1,0000 

.0000 

0.0000 


3.62 

1.12 


-.0465 

,9666 

,0862 

-.0063 

. . . 

3.66 

1.13 


-.0492 

.9610 

.0955 

-.0074 

. , 

3.85 

1.30 

. ... 

-.0735 

,8155 

.2895 

-.0315 


3,56 

1.30 


-.0735 

.8165 

.2895 

-.0316 


4.04 

1.43 


-.0699 

,6570 

.4655 

-.0527 


3,88 

1.49 



-.0637 

.5762 

.5487 

-.0612 


3.60 

1.58 


-.0612 

.4517 

.6701 

-.0706 


3.98 

1.75 


-.0234 

.2266 

.8672 

-.0703 

. . , 

3,70 

1.87 


-.0074 

.0956 

,9610 

-.0492 


3.87 

1.03 


.0000 

-.0141 

.9978 

.0168 

-.0004 

3,58 

1.22 



-.0669 

.8960 

.1908 

-.0189 

3,69 

1.30 

, 


-,0735 

.8155 

.2895 

-.0315 

3.29 

1.33 



-.0741 

.7417 

.3289 

-.0365 

3,56 

1.43 



-,0699 

.6570 

.4655 

-.0527 

3.27 

1.55 



-.0557 

.4933 

.6304 

-.0681 

3.10 

1.57 

. . , 


-,0527 

.4655 

.6570 

-.0699 

3.23 

1.78 



-.0189 

.1908 

,8950 

-.0669 

2.74 

1.80 



-.0160 

.1680 

.9126 

-.0640 

2.98 

1.93 


* t % 1 

-.0023 

.0443 

,9883 

-.0303 

2.55 


^1 

^2 

^3 

Z4 

-^5 

^6 

y 


Conventional regression symbols 


1 

From table A-1, 



^ ^ , values 03 coetlicients ' 

Base-point Scale for base points 

positions variate 

[nx 0 YBl YB2 YB^ YB4 YB5 7^6 YB7 YBS K 

.6 

.8 ... 


2.06 

0.0048 

- 0.0373 

0.9949 

0.0429 

- 0.0054 

0.0001 



1.50 

2.08 

.0062 

-.0483 

.9907 

.0588 

-.0075 

.0001 



1.69 

2.29 

.0149 

•-.1144 

.8566 

,2842 

-.0454 

.0040 



2.02 

2.25 

.0143 

-.1090 

.8948 

2332 

-.0360 

.0028 



2.20 

2.32 

.0151 

-.1163 

.8246 

.3246 

-.0628 

.0050 



2.50 

2.50 

,0117 

-.0977 

,5859 

.5859 

-.0977 

.0117 



2.50 

2.50 

.0117 

-.0977 

.5859 

.5859 

-.0977 

.0117 



2.70 

2.57 

,0092 

-.0813 

.4826 

.6858 

-.1101 

.0138 



2.89 

2.83 

.0010 

-.0202 

.1421 

.9536 

-.0882 

,0115 



3.25 

2.87 

.0005 

-.0139 

.1026 

.9739 

-.0726 

.0095 



3.58 

3.00 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 



3.62 

2.12 


.0089 

-.0681 

.9780 

.0933 

-.0125 

0.0004 


3.66 

2.13 


.0095 

-.0726 

.9739 

.1026 

-.0139 

,0005 


3.85 

2.30 


.0150 

-.1152 

.8462 

.2975 

-.0478 

.0043 


3.55 

2.30 


.0150 

-.1152 

.8462 

.2975 

-.0478 

,0043 


4.04 

2,43 


.0138 

-.1101 

.6858 

.4826 

-.0813 

.0092 


3.88 

2.49 


.0121 

-.0997 

.6005 

.6713 

-.0955 

.0114 


3.60 

2.58 


.0088 

-.0787 

.4678 

.6995 

-.1114 

.0140 


3.98 

2.76 


,0028 

-.0360 

.2332 

.8948 

-.1090 

.0143 

. 

3.70 

2.87 


.0005 

-.0139 

.1026 

.9739 

-.0726 

.0095 


3,87 

2.03 



.0026 

-.0194 

.9988 

.0207 

-.0026 

0.0000 

3.58 

2.22 



.0135 

-.1028 

.9197 

,1972 

-.0296 

.0020 

3.69 

2,30 



.0150 

-.1162 

.8462 

.2975 

-.0478 

.0043 

3.29 

2.33 



,0151 

-.1166 

.8132 

.3383 

-.0554 

.0054 

3.66 

2,43 



,0138 

-.1101 

.6858 

.4826 

-.0813 

,0092 

3.27 

2.55 



.0099 

-.0862 

.5122 

,6579 

-.1071 

.0133 

3.10 

2,67 



.0092 

-.0813 

.4826 

.6858 

-.1101 

.0138 

3.23 

2.78 



.0020 

-.0296 

.1972 

.9197 

-.1028 

.0135 

2.74 

2.80 



,0016 

-.0256 

.1744 

.9344 

-.0976 

,0128 

2.98 

2.93 



,0001 

-.0064 

.0508 

,9930 

-.0429 

.0055 

2.55 


Xi 

^2 

^3 

^4 

^5 

^6 

Xj 

^8 

Y 


Conventional regression symbols 


Vrom table A * 3 . 


NOTES FOR PROGRAMS B-1 AND B-2 

AA(I) 

Matrix to be inverted. 

X(I} 

Input conventional re- 

MINVIAA.N.D.B.M) 

System library routine 
for matrix inversion. 


gression data from 

BID 

Regression coefficients 

SSX(I) 

tables B-2 and B-3. 
Sums of squares and 

INTFOR{V,Y,YA.MS) 

or base points. 

See program A-1. 


products. 

INTSIX(V,Y,YP,MS) 

See program A-3, 




PROGRAM B-1.- — Four-point smoothing 


dimension X(30,10) ,SSX(10,10) iL( 10),M(10) »YINT(60) , YAINTtfeO) 
1,AA(6,6) 

COMMON/XX/ TAB(100» 

READ <5,10 00) ((X (I, J> ,J=1,7) ,1=1,30) 

1000 F0RMAT{7F5.4) 

00 1002 1 =1,10 
DO 1002 J=l,10 
1 002 SSX(I,J) = 0.0 
00 1003 1=1,7 
00 1004 J=I,7 
no 10 0 5 K = l,30 

1 0 05 SSX (I,J) = SSX(I,J)+X(K,I)<‘X(K,J) 

1004 CONTINUE 
1003 CONTINUE 

DO 1006 1=1,7 
DO 1006 J = I,7 
10 06 SSX (J ,I) = SSX(I, J) 

NRITE (6,1010) ( (SSX (I, J) ,1=1,7) ,J=1, 7) 

1010 FORMAT(* *,(7F12.8)) 

DO lOZO 1=1,6 

DO 1020 J=l,6 
1020 AA<I,J)=SSX(I, J) 

CALL MINV(AA,6,0,L,M) 

WRITE <6, 10 11) ((AA (I,J) ,1 = 1,6) ,J=1,6) 

1011 FORHAT(* »,(6F13.8)) 

DO 1007 1=1,6 

8 (I ) = 0 . 0 
DO 1008 J=l,6 

1008 8 (I)=B(I) +AA(I,J)*SSX(J,7) 

1 0 07 CONTINUE 

WRITE(6, 1011) (B(I) ,1 = 1,6) • 

TAR(1)=0.4 
TAB(2)=0. 2 
DO 1012 1=3,8 

1012 TAB(I)=B(I-2) 

V = 0. 6 

MS=-1 

DO 1013 1=1,31 

CALL INTFOR(V,Y,YA,MS) 

YINT(I)=Y 
YAINT(I)=YA 
V=V+0, 02 

1013 CONTINUE 
YAINT(1)=0.0 

WRITE(6,1014) (I, YINTd) ,1 = 1,31) 

1014 FORNAK* INTERPOLATED FUNCTION*/ (8 ( 15, FIO . 6) ) ) 

WRITE (6, 10 15) (I,YAINT (I) ,1=1,31) 

1015 FORMAK* INTERPOLATED INTEGRAL*/ (8 (15, FIO . 6) ) ) 

STOP 

END 



PROGRAM B-2. — Six-point smoothing 


OIMENSION X(30',10) , SSX ( 10 , 10)j L'( lOV, M( 1 0) ,B(iO) ,YINf(60) ,YPINT{60) 

liAAte,*)) 

COMMON/XX/ TAR(IOO) 

^EAO(5,10 00)<(X(I, J) ,J=:1.,9> ,1 = 1 *30) 

ITO'O FORMAT (gp 5. A) 

GO 1002 1=1,10 

■■ 00" 1002 J = l,10 ■ ■■ 

1002 SSX(I,J)=0.0 
00 1003 1=1,9 
00 lOOA J=I,9 

90 1005 K = 1 ,30 ■ ■ 

1005 SSXd, J) = SSX(I,J)+X(K,I)»X(K,J) 
lOOA CONTINUE 

1003 CONTINUE 

no 1006 1=1,9 ■ ' 

no 1006 J=I,9 

1006 SSX (J,I>=SSX(I, J) 

WRITE (6, 1010) ((SSXd, J),I = 1,9) ,J=1,9) 

1010 F0RMAT(* ♦,(9F12.8)) 

00 1020 1=1,8 

00 1020 J = l,8 ' •- 

1020 AAd,J) =SSXd, J) 

CALL MINV(AA,8,0,L,M) ' 

WRITE (6, 1011) ((AA(I,J) ,1 = 1,8) ,J = 1,8) 

1011 FORMAT(^ •,(8F13.8)) ’ ' 

no 1007 1=1,8 

nd)=o.o ' - 

00 1008 J=l,8 

1008 9 (I)=Bd) +AAd,J)^fSSX(J,9) ' " ‘ 

1007 CONTINUE 

WRITE(6, 1011) (Bd) ,1=1,9) - 

TAB(1)=0. A 
TAB(2) = 0,2 
DO 1012 1=3,10 

1012 TABd) = B(I-2) 

V=0,8 

MS=-1 — 

00 1013 1=1,31 

CALL INTSIX(V,Y,YP,MS) ' 

YINTd)=Y 

YPINTd) = YP - 

V=V+O,02 

1013 CONTINUE 


WRITE(6,1014) d,YINT d) ,1 = 1,31) 
lOU FORMATd INTERPOLATED FUNCTION V (8(15, KlO. 6)) y 
WRITE(6, 1015) d,YPINT(I) ,1=1,31) 

^ INTERPOLATED DERIVATIVE*/ (8{l5,>ro , 6)' )) 


END 
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